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We present a general theory for semiconductor polarons in the framework of the Frohlich in- 
teraction between electrons and phonons. The latter is investigated using non-commuting phonon 
creation/annihilation operators associated with a natural set of non-orthogonal modes. This setting 
proves effective for mathematical simplification and physical interpretation and reveals a nested 
coupling structure of the Frohlich interaction. The theory is non-pert urbative and well adapted for 
strong electron-phonon coupling, such as found in quantum dot (QD) structures. For those particu- 
lar structures we introduce a minimal model that allows the computation and qualitative prediction 
of the spectrum and geometry of polarons. The model uses a generic non-orthogonal polaron basis, 
baptized "the natural basis" . Accidental and symmetry-related electronic degeneracies are studied 
in detail and are shown to generate unent angled zero-shift polarons, which we consistently elim- 
inate. As a practical example, these developments are applied to realistic pyramidal GaAs QDs. 
The energy spectrum and the 3D-geometry of polarons are computed and analyzed, and prove that 
realistic pyramidal QDs clearly fall in the regime of strong coupling. Further investigation reveals 
an unexpected substructure of "weakly coupled strong coupling regimes", a concept originating 
from overlap considerations. Using Bennett's entanglement measure, we finally propose a heuristic 
quantification of the coupUng strength in QDs. 



I. INTRODUCTION 

Quantum structures (QSs), such as quantum dots 
(QDs), are sophisticated solid-state pieces, vital for fun- 
damental research and novel applications in quantum op- 
tics and quantum informatics. Today, QDs find techno- 
logical use in QD lasers^ infrared photodetectors^ sin- 
gle photon sources^ or markers in biology. Cutting- 
edge research features QDs as medical fluorophores for in 
vivo detection of cell structures such as tumors.^ Other 
promising applications for QDs are solar cells^ and opti- 
cal telecommunication.^ The most exciting, yet challeng- 
ing expectation relies in the use of QDs as qubit holders 
and gates for quantum computation.^ For fundamental 
science, QDs are among the few systems allowing con- 
trolled experiments with single energy quanta giving di- 
rect access to controlled quantum entanglement and cor- 
relations. 

Due to their extreme carrier sensitivity, much interest 
in QDs relates to carrier relaxation and excitation pro- 
cesses mediated by various interactions, such as carrier- 
carrier, carrier-photon and carrier-phonon interactions. 
As for carrier-phonon interactions, early pert urbative ap- 
proaches with acoustic phonons resulted in the bottleneck 
concepti^^^^ Theses perturbative results predict in- 
efficient carrier relaxation for a large class of small QDs. 
Although experimentally verified in certain casesi^^I^, 
these predictions failed in many other tests ^i^iii A defi- 
nite progress came with non-perturbative investigations 
of the deformation potential and Frohlich interaction, re- 
vealing the existence of a strong coupling regime, which 
is out of reach of perturbative approaches and allows 
efficient carrier relaxation through acoustic and optical 
phonon dynamics respectively^i^ii^i^Si^ This led to the 
new concept of quantum dot polarons (QDPs), which 



are non-separable fundamental excitations determined 
by the carrier-phonon interaction. Within the approx- 
imation of monochromatic LO-modes for the Frohlich 
interaction, electrons only couple to a finite number of 
lattice modes as analytically explained through an alge- 
braic decomposition introduced by Stauber et al?'^ Their 
procedure constructs an orthonormalized basis of rele- 
vant lattice modes from the finite set of phonon cre- 
ation/annihilation operators naturally appearing in the 
Frohlich Hamiltonian. This leads to a numerically solv- 
able model of QDPs^^, which can be viewed as an exten- 
sion of the work by Ferreira et al?^ 

In this work, the polaron problem is tackled from a 
different viewpoint: the full electron-phonon Hamilto- 
nian is reformulated in terms of non-orthogonal modes, 
which naturally span all coupled and uncoupled crystal 
vibrations. The non-orthogonal structure is preserved 
from the beginning to the end and exhibits undisputable 
advantages for computation and physical understanding. 
General analytical results applicable to any type of semi- 
conductor QS are derived in this framework. They are 
subsequently applied to peculiar pyramidal Csy GaAs 
QDs, but the same theoretical scheme could be applied 
to any other semiconductor QD structure, e.g. zincblende 
InAs QDs with C2v symmetr}^ or Wurzite GaN QDs 
with high Csy or C^y symmetry. 

Section HIl considers a general QS populated by an ar- 
bitrary number of bound electrons and phonons. We 
first introduce a set of non-orthogonal LO-modes, which 
spans all the LO-modes appearing in the Frohlich inter- 
action. From there we derive two decoupled subalgebras 
of non-commuting phonon creation/annihilation opera- 
tors, which separate the quantum structure in a subsys- 
tem of bound polarons and a subsystem of uncoupled 
modes pi Bp . The theory culminates in a non-trivial 
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nested coupling structure of the Frohlich interaction, 
which has important consequences when working with 
any finite number of phonons (|II Cp . In section HTH we in- 
troduce a minimal non-perturbative model (one electron, 
one phonon) particularly suitable for the crucial case of 
QDPs. We provide an explicit non-orthogonal polaron 
basis, baptized the "natural basis" pil A\i . It provides 
a detailed interpretation of the geometries and spectra 
of low-energy QDPs. We also investigate additional sim- 
plifications resulting from electronic degeneracies piIBp 
and group theoretical considerations of dot symmetries 
(jlll Cp . The theory concludes with some key aspects of 
the three-dimensional (3D) numerical code pll Dp . which 
comprises an adaptive irregular space discretization for 
computing the Frohlich matrix elements. 

In a second part, the minimal model is applied to 
realistic pyramidal GaAs QDs with Csy symmetryj^i 
Section IIVI presents the 3D-geometries of the QDPs 
and their spectrum, throughout using group theory 
(|IV A( IIVB| lIVEp . Explicit comparison with perturba- 
tion theory confirms the existence of a strong coupling 
regime. Surprisingly, we find significant numerical evi- 
dence for a peculiar substructure inside the strong cou- 
pling regime. This leads to the concept of "weakly cou- 
pled strong coupling regimes" pV Cp . which can be un- 
derstood in terms of overlap between confined electrons 
and coupled modes. Using Bennett's entanglement mea- 
sure, we further present a useful alternative character- 
ization of the strong electron-phonon coupling in QDs 
pVDl) . 

In section lYl we report on the general polaron proper- 
ties that could be expected in other QD systems. Section 
|Vl] concludes the theory with a short review, and helpful 
derivations are provided in the appendix I VII[ 



Hamiltonian iJ*^*, 

(1) 

(Unity operators and tensor products have been omit- 
ted.) a^cr, cijj^a fermionic annihilation and creation 
operators of confined conduction electrons, with /i label- 
ing an orthogonal set of stationary wave functions and 
a being the spin index. The scalars are the free elec- 
tronic energies, which are independent of a in the ab- 
sence of magnetic fields. 6q, 6^ are the bosonic annihi- 
lation and creation operators of phonons associated with 
the LO-plane waves ^q(x) = Y^2/Fe*^'^, where V is the 
quantization volume, clq = ^LO is the phonon en- 
ergy assumed independent of q (monochromaticity), and 
Mjj^jj^'d are the Frohlich matrix elements — 

^""'^ ^\]^^ {'^'t:^^ (2) 

where Sstat and Soo are the static and high frequency di- 
electric constants and '0/i(x) are the (one-particle) elec- 
tronic wave functions. Since the Hamiltonian ([T]) is de- 
coupled and symmetrical in spin degrees of freedom, we 
shall from here on omit the spin indices a. 



B. Subsystem of Quantum Structure Polarons 



II. NON-ORTHOGONAL THEORY FOR 
POLARON STATES 

In this section, we present a theory for polar semicon- 
ductor QSs, e.g. dots, wires or wells, in which the carrier 
evolution is reasonably described by Frohlich interactions 
with monochromatic LO-modes. The QSs can contain an 
arbitrary number of electrons (within the limitations in- 
duced by the Pauli exclusion principle) and an arbitrary 
number of phonons. The conservation laws exhibited by 
the interaction Hamiltonian allow straightforward gener- 
alizations to exciton-polarons or even polarons associated 
with bigger electron/hole complexes. 



A. Polaron Hamiltonian in Quantum Structures 

The model's evolution is dictated by a Hamiltonian 
composed of a free evolution term and the Frohlich 



We shall now apply a non-orthogonal linear transfor- 
mation to the operator basis {^q, ^q} in order to reveal 
two decoupled physical subsystems, the subsystem of 
"Quantum Structure Polarons" (QSPs) and the subsys- 
tem of "Uncoupled Phonons" (Uphs). This conceptual 
separation will be refiected in a tensor product decom- 
position of the representative Hilbert space. 

The matrix elements ([2]) can be considered as discrete 
3-dimensional functions of q. They obey relations of lin- 
ear dependance, as can be seen by choosing the electronic 
wave functions ^l^fj,{x) real (always possible), in which 
case M^^/q = M^/^q. If there are N orthogonal electron 
states /i, the number of such relations is A^(A^ — 1)/2. The 
remaining N{N -\- 1)/2 matrix elements show no obvious 
relations of linear dependance, and we shall temporarily 
assume that there are the only N{N — l)/2 independent 
relations of linear dependance. The theory remains valid 
in the case of additional linear dependencies such as dis- 
cussed towards the end of this subsection. 

The structure of the Frohlich interaction implies that 
the number of linearly independent matrix elements 
^/i/i'q equals the number of linearly independent lattice 
modes that appear in the interaction term. This can be 
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seen explicitly, when reformulating the interaction as 

= E + h.c. (3) 

^MM'^ 7^^M^^,q6q = ^L^^,q6q (4) 



where \Jfifi'\'^ = Xlq I^MM'qP quantizes the electron- 
phonon coupling strength, with J^^/ chosen as a posi- 
tive real. The relations of linear dependance among the 
matrix elements M^^/q trivially translate to B^^^, — Sj^,^ 
and B^^' = B^>^. The remaining N{N ^1) /2 linearly in- 
dependent phonon operators shall be scanned by a unique 
pair index A = {/i, ji'} = {/i^ fi}. 

The operators {5a, B\} annihilate and create "coupled 
phonons" , that is quanta in terms of a harmonic oscillator 
in modes susceptible to interact with electrons via the 
Frohlich potential. Using (|4]), the wave functions of those 
modes are given by the inverse Fourier transforms 



q 



iq x r * 
6 ^Aq 



(5) 



which are manifestly localized in the quantum structure. 
The set of all modes Sa(x) is non-orthogonal as empha- 
sized by the non-diagonal scalar product matrix and the 
non-diagonal commutator of the corresponding operators 
{Bx.bI}. Both follow directly from © and 



A 



Bx,Bl 



Aq 



(6) 



(Round brackets represent the scalar product relative 
to the quantization volume V.) For reasons of physi- 
cal interpretation and mathematical simplicity, we skip 
a possible orthonormalization and preserve the non- 
orthogonality for the rest of the theory. 

In order to express the full Hamiltonian in terms of 
the new operators {5|,5a}, we need to complete them 
by an operator-set {B^, 5q} generating the orthogonal 
complement {Sq(x)} of the coupled modes {Sa(x)}. We 
choose a linear transformation. 
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'qq'^q 



= 5Zcqq/^q(x) (7) 



A natural and sufficient condition for the coefficients Cqq/ 
writes ^q|0) = (li^/j, — 7^) 6^|0), where |0) is the phonon 
vacuum, ti ph the unity on the subspace of one phonon, 
and V is the orthogonal projector on the sub-subspace of 
coupled one-phonon states 'uectj^J I O)}. Thus (lip/i — 7^) 
projects on the one-phonon sub-subspace of uncoupled 
modes, and {^q} creates quanta accordingly called "un- 
coupled phonons". An explicit derivation of the coef- 
ficients Cqq/ is provided in appendix IVII A[ From this 
explicit form it follows that the modes {Sq(x)} are 
also mutually non-orthogonal, which again translates to 



a non-diagonal commutator of the corresponding cre- 
ation/annihilation operators {5^, 5q}, 



^q' ~ ('^q'^q') ~ ^^^qq^^q'q"' 



(8) 



Indeed, the modes {Sq(x)} constitute an overcomplete 
set, according to the N{N + l)/2 relations of linear de- 
pendence 



(9) 



However, it is important to note that all coupled modes 
{Sa(x)} are orthogonal to all uncoupled ones {Sq(x)}, 
as emphasized by the following commutators and scalar 
products 



[Bx,Bl] = =0 



(10) 



The transformations (|4]) and ([7j) constitute a non- 
orthogonal mapping {6^} The inversion 
{^q} ^ {^A' ^q} unique due to the over complete- 
ness of {5^}. A suitable form, consistent with (|4]) and 
([7j), is given by 



^q = 4 + EWA"')AA'5l 



(11) 



AA' 



This allows us to express the phonon number operator in 
terms of the new operators. 



E^^^. = E(^"')av51^v + E^^^. 

q AA' q 

Finally, the full Hamiltonian ([T]) transforms to 
H° ^ ^e^ata^ + eioE(A"')AA'^A^^' 



AA' 



■^Uph 



(12) 

(13) 
(14) 
(15) 



(unity operators and tensor products have been omitted). 
The fundamental commutators (p!Q|) imply the commuta- 
tor 



(16) 



The latter defines a unique separation in two physical 
subsystems, expressed by the tensor product decomposi- 
tion 



n = H^^^ (g) H^P^ 



(17) 



such that H^^^ acts trivially in H^p^ and H^p^ acts 
trivially in H^^^. [1^ assumes the bosonic sym- 
metrization of the phonon subsystem.] The subsystem 
represented in 7i^^^ consists of electrons and coupled 
phonons associated with a finite number of A^(A^ + l)/2 
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linearly independent modes {Sa(x)}. The stationary 
states (i.e. eigenstates of H^^^) are likely entangled in 
electronic and phononic coordinates and will be referred 
as to "quantum structure polarons" (QSPs). In contrast, 
the subsystem of "uncoupled phonons", represented in 
^Uph^ is a pure phonon-system associated with infinitely 
many uncoupled bulk modes {Sq(x)}. Each such mode 
evolves trivially under the phonon number operator, and 
thus the quantum structure problem drastically reduces 
to solving H^^^ inside . 

This theory remains valid if the Frohlich matrix ele- 
ments exhibit other linear dependencies than M^^t^^ = 
^/i'/iq (and their linear combinations). Indeed, if there 
are N' linearly independent matrix elements, N' < 
N{N -\- l)/2^ it suffices to redefine the index A such as 
to label only the corresponding independent operators 
B^^^'. The derivation above stays valid with this redefi- 
nition, if any number N{N ^1) /2 is replaced by N' . (e.g. 
the number of linear dependencies among the uncoupled 
modes Bq will be reduced to N' ^ etc.) 

It is worth noting that the Hamiltonian H^^^ mani- 
festly conserves the number of electrons, i.e. 







(18) 



This conservation law implies the existence of one cou- 
pled mode that only couples to the electron number op- 
erator, such as shown by Stauber et al.^^ In contrast to 
their choice, we decide to keep this particular mode in 
the system of QSPs. Indeed, even though this mode does 
not affect the overall electron dynamics, it is located in 
the quantum structure and evolves through the creation 
an reannihilation of intermediate electrons. Therefore, 
its stationary solutions are Glauber-coherent states, very 
different from the stationary phonon-number states of 
uncoupled modes. 



C. Nested Coupling Structure 



reveals that coupling structure can be identified individu- 
ally for each fixed number of electrons without loss of gen- 
erality. For the rest of this section, we shall thus restrict 
our considerations to some fixed number of electrons m 
(m > 1), and take the subspace VP^^ as restricted to m 
electrons. Second, we note that H^^^ does not couple 
orthogonal spin states, and hence the coupling structure 
can be investigated with all electrons having the same 
fixed spin a. As H^^^ acts identically on all values of a 
spins can be generally neglected (as in the previous sec- 
tion). Third, we use the property that the Frohlich op- 
erator H^^^ affects phonon numbers by one unit. Hence, 
it is useful to decompose VP^^ in subspaces associated 
with different numbers of phonons /c. 



QSP 



k=0 



n 



QSP 



vect 



{at. 



■B^ ,\0)] 



(20) 



The index list {i^^, /i^, /i^} goes over all combinations of 
electronic indices, such that Vi ^ \l i ^ j (Pauli ex- 
clusion principle). Here, |0) = \Oelectrons) \Ophonons) 

denotes the polaron vacuum. 

According to the coupling rule ([19]), the subspace cou- 
pled to H^^^ is given by 

vect {e-'^^""*|^) : Vt, 1^^) G Tif} (21) 

To pinpoint a particularity in the coupling between 
H^^^ and its "inferior neighbor" H^^^ ^ we shall tem- 
porarily restrict the phonon Fock space to at most p 
phonons {p > 0). This implicitly requires a trunca- 
tion of the Hamiltonian H^^^ equivalent to imposing 
B{\ip) = \ip) e n^^^. We define H^^^ as the 
sub-subspace of TL^^^ coupled to Tipfi within this re- 
striction. Departing from ([2T]) with k = p — 1, H^^^ can 
be simplified to (derivation in appendix IVII Bp 



In this section, we shall uncover a nested structure in 
the Frohlich coupling. This structure implies in partic- 
ular that certain states differing by one phonon (e.g. a 
state with one phonon and a state with two phonons) 
are exclusively coupled via intermediate higher order 
states (e.g. a state with three phonons). This non-trivial 
coupling structure provides some intuition for the form 
of stationary states and implies a rule to truncate the 
Hilbert space if the polaron problem is restricted to a 
finite number of phonons. 

In the following non-pert urbative analysis, two states 
\ (p) are called "coupled" if the evolution of one state 
develops a non-vanishing projection on the other, i.e. 
(^ilj\exp{—iH^^^t)\(p) 7^ for at least one t. Thus the 
subspace Ha coupled to a subspace Hb is given by 



(19) 



In order to identify a coupling structure, we first use 
the conservation of the number of electrons [eq. (p!8|)]. It 



vect < e 



(22) 



where is the free evolution ([T4j) and 

Xl/i/i' ^MM'^iL^M'^/i/i' denotes the phonon creating, pctxu w± 
the Frohlich interaction (|3]). In physical terms, ([22]) ex- 
presses that an electron-phonon state | (/:)), initially con- 
taining p — 1 phonons, evolves towards a superposition 
involving a certain p-phonon state (by Frohlich interac- 
tion). The latter is generally not an eigenstate of and 
its free evolution can span a whole p-phonon subspace 
coupled to the initial state \(p). For further simplifica- 
tion we decompose H!^^\(p) in eigenstates of i^^. 
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(23) 



where 7 labels the eigenspaces of H'^ inside H^^^ ^ and 
are the orthogonal projectors on all these eigenspaces. 
As projects on a p-phonon subspace and \(p) is 
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a (p — l)-phonon state, we can safely replace V^H!^^ On the other hand, there may be stationary polarons 
by T^^iJ*^*, for annihilates the phonon annihilat- 



ing part of the interaction H^^^ 



in (|22|) gives 



Invoking the relation 
(f) and substituting 



TY^^^ = vect i^V^H'^'l^) : V7, \^) e (24) 

(since e~*^^^/^ for different 7 are linearly independent 
functions of t.) 

The eigenspace projectors act trivially on the sub- 
system of lattice modes populated by p phonons, since 
all j9-phonon states are degenerate (monochromaticity as- 
sumption). As for the electron subsystem (here consid- 
ered as non-degenerate), the different eigenspaces can be 
labeled as 7 = {i^i, . . . , i/^}, where {. . .} denotes an un- 
ordered set and vi i 3 - (Spin indices were omit- 
ted according to the introduction of this section.) The 
electronic part of the projectors V^^ can then be expressed 
as 



' ' ' l^electrons electrons \ ^Vrr 



(25) 



Substituting this expression in (|24|) , allows to express 
H^^^ with explicit basis vectors. After rearrangement 
and substitution of indices, we find 



vectl a 



t 



•••B 



t 







)} (26) 



where /ii = i^i, and 7^ Ui Vi = 2, . . . , m. Expression 
([26]) shows that H^^^ is necessarily a subspace of H^^^ . 

We conclude that if the number of phonons is limited 
to p (p > 1), the subspace Hpfi couples to H^^^ , but 

not to its orthogonal complement H^^^^ C H^^^ . 

In conclusion, if for physical or computational reasons 
the model is truncated to a finite number of phonons p, 
then the subspace H^^^ must be restricted to H^^^ . 
Otherwise non-physical polarons would appear (con- 
tained in H^^^^), that would seem uncoupled and thus 
unshifted relative to the free spectrum. Such a precau- 
tion was apparently not taken in previous works. The 
particular truncation H^^^ H^^^ also represents an 
analytical and computational simplification. 

If we release the temporary assumption of a finite 
phonon number p (or if we take k < p)^ the following 
statement holds: /c-phonon states in H^^^ ^ do not di- 
rectly couple to {k — l)-phonon states, but can only cou- 
ple to {k — l)-phonon states via intermediate {k + 1)- 
phonon (and higher order) states! In a perturbative ap- 
proach, these particular couplings would first appear in 
the third order of the interaction term. Direct couplings, 
i.e. couplings that do not involve intermediate higher- 
order states, are represented by the arrows in Fig. [H 
This nested structure provides some insight in the form 
of stationary polarons (which generally superpose states 
with different phonon numbers). For example, stationary 
superpositions of states from H^^^ ^ and H^^^ neces- 
sarily involve a strong contribution of states from H^^i • 



made of states from H^^^ and H^^^ with only a minor 
contribution of states from H^^^ . 




etc. 



phonons 1 phonon 2 phonons 



3 phonons 



FIG. 1: Nested coupling structure exhibited by the Frohlich 
interaction. Arrows represent couplings that do not involve 
intermediate states with a higher number of phonons (see 
text). Solid circles contain the subspaces Ti^^^ for differ- 
ent /c's. Dashed circles enclose the subspaces H^^^ , whereas 
hatched zones show their orthogonal complements Ji^^^ ^ 
within Ti^^^ . The solid gray filling indicates the subspace 
7-^* spanned by the natural basis (|28)) . 



III. ONE-ELECTRON/ONE-PHONON MODEL 
OF QDPS 



In the framework of the general non-orthogonal the- 
ory developed above, we shall now propose a minimal 
non-perturbative model for polaron states in quantum 
dots (QDs). The general quantum structure considered 
so far, is now specified as a quantum dot: QS^QD 
and QSP^QDP. In such zero-dimensional systems, the 
monochromaticity assumption, crucial for the present 
theory, is fairly precise for the relevant modes (i.e. wave- 
lengths comparable to the dot size and thus long com- 
pared to the atomic spacing). The model assumes a sin- 
gle electron (m = 1) populating different levels while 
coupling to at most one phonon {p = 1). It allows to ap- 
proximate the shifts of the first polaron levels, which are 
typically populated at low temperatures, although there 
may be additional effects arising from acoustic phonons 
like dephasing effects. 

In the next three subsections, we subsequently inves- 
tigate QDs with non-degenerate electron levels pil Ap . 
with accidental degeneracies pil Bp . and with symmetry- 
related degeneracies pil Cp . For each case, we develop 
a simple non-orthogonal polaron-basis S*, baptized the 
"natural basis", which spans the relevant Hilbert space 
7Y*. A similar formalism could be developed for holes 
(although acoustic phonons may have to be taken into 
account there), or for any many-particle complex such 
as an exciton-, a trion- or a biexciton-based quantum 
dot polaron. Stauber and Zimmermanr^ showed that a 
correction term must be introduced in the case of non- 
neutral complexes. 
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A. Natural Basis 

We first consider a non-symmetric QD witii N non- 
degenerate electron levels ja. Accordingly there are 
A^(A^+l)/2 linearly independent coupled modes, spanned 
by the operators sj^^^,. By virtue of the coupling struc- 
ture developed in section III C[ the coupled regime of the 
one-phonon model is properly represented by the sub- 
space 7Y* corresponding to gray filling in Fig. [TJ It writes 

= n^^^ © Hf^^ (27) 

A vector set S*, such that 7i* = vectB* is directly ob- 
tained from ([26|) . 

r = {|M;0),B^^,|/i;0) = ^M;^,q|/i;q) V/i,/i^} (28) 

q 

where we used the short hands |/i;0) = a]^|0) and 

|/i;q) = a^bl^lO), with |0) = \Oelectrons) © \Ophonons) 

being the polaron vacuum. The vectors in ([28|) are gener- 
ally non-orthogonal but linearly independent and S* will 
be called the "natural basis". All natural basis states 
are eigenstates of H^. For each electronic level ja there 
is one natural basis state with zero phonons (free energy 
e^) and there are N natural basis states with one phonon 
(free energy + eLo)- Since there are N electronic levels 
/i, the dimension of the relevant subspace 7Y* writes 

dim{n'^) = card{B*) = N{N + 1) (29) 

The requirement to reduce the one-phonon subspace 
^QDP H^^^ (section III Cp reveals the simplifying 
feature that many product states of electron states and 
coupled phonons are irrelevant for the polaron structure 
(e.g. ^i^/xsl/^i^O) i ^*)- Therefore, the number of QDPs 
only scales as N'^ and not as (!), which one might ex- 
pect from the number N of dot electron states and the 
number oc N'^ of coupled modes. 

Fig. [2] shows the qualitative QDP spectrum in the case 
of a QD with only three non- degenerate electronic states. 
Gray bars denote additional QDPs that would appear in 
an extended model including the interaction with two- 
phonon states. Those are associated with free states in 
^QDP± j^.^ ^j^^ orthogonal complement of H^^^ inside 
^QDPy The connections between free levels and QDP 
levels (Fig. [2]) are an important outcome of the natural 
basis. They indicate the free levels from which specific 
QDPs arise, if one could gradually introduce the Frohlich 
interaction. This picture allows a prediction of spectral 
changes under dot size variation, since one can generally 
assume that shifts increase when approaching a resonance 
of the Frohlich interaction {i.e. Ae = e^o)- 

In conclusion, a complete non-orthogonal polaron basis 
called the "natural basis" ([28|) has been introduced. It 
provides a mean for physical understanding of polaron 
spectra involving low phonon numbers, and constitutes a 
simplifying and powerful computational basis (see IIII D[ 

IIYI). 




FIG. 2: Qualitative structure of the polaron spectrum in the 
case of three bound, non-degenerate electronic levels. Each 
gray bar indicates 3 additional polaron levels that would result 
from interactions with two-phonon states. 



B. Accidental Electronic Degeneracies 

This section and the next one point out additional 
simplifications in the case of electronic degeneracies. 
In particular, if a non-degenerate electronic spectrum 
(e.g. Fig. [2]) becomes partially degenerate, for example by 
specific dot size adjustment, not only certain QDPs may 
become degenerate, but some of them will analytically 
align with free levels. We shall call such states "zero-shift 
polarons" and show that they are nothing but uncoupled 
states, susceptible to become QDPs as soon as the degen- 
eracies are lifted. Thus the relevant Hilbert space 7Y* can 
be further reduced, such that spurious zero-shift polarons 
are automatically eliminated. 

In order to label accidental degeneracies, the electron 
index is now expressed as /i = (r, z), where r = 1 . . . n < 
N is an energy index and i = 1 , . . . , a degeneracy 
index. The eigenspaces of inside the one-phonon 
subspace Tif^^ are indexed by 7 = r, and the or- 
thogonal projectors = Vr on these eigenspaces write 
Vr = X^f^ \r^i; q)(^r^i; q\. Substituting these projec- 
tors in ([24|) allows to write the relevant subspace 7Y* as 

(30) 

with H^^^ = vect{VrH'''^\^]{)) Vr,/i} 

Expressing Hq^^ , Vr^ and in basis vectors |/i;0) 

and q), naturally provides a basis of 

= {|r,z;0),f^ ^M^,,,,.^|r^^^q) Vt,/,z} (31) 

i' = l q 

Its dimension is 

dim{n*) = card{B*) = N{n + 1) (32) 

where N is the total number of orthogonal electronic 
states in the dot and n < is the number of distinct 
electronic energies (n=N would be the non-degenerate 
case.) We note, that even though the number of polarons 



7 



is smaller in the degenerate case, the number of modes in- 
volved remains the same. Only the number of accessible 
product states [electron) |phonon) is reduced. This can 
be seen from (jSTJ), which, for a degenerate level r, yields 
entangled states similar to YlP' i){T' '^^)^ 

As in the previous section, the natural basis (|3T]) pro- 
vides a qualitative prediction of the polaronic spectrum 
and associates each polaron level with a free level (see 
example Fig. [3]). In particular, we emphasize that the 
highest free level in the figure only yields 3 orthogonal po- 
laron states and not 6 as one might expect from pulling 
together the two uppermost free levels in Fig. [2l The 
particular case of symmetry related degeneracies is now 
addressed in the next subsection. 



QD potential 
and electronic 
wave functions 




Free 
levels 



QDPs 



2x 



1 X 



1 X 
1 X 
1 X 
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1 X 
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FIG. 3: Qualitative structure of the polaron spectrum in the 
case of two electronic levels, one of which is twice degenerate. 
Each gray bar indicates 3 more polarons that would result 
from interactions with two-phonon states. 



C. Symmetrical Quantum Dots 

Additional degeneracies and simplifications may be ob- 
tained in the case of QDs invariant under a set of symme- 
try operations, generally described by the group of such 
operations i. e. [H^g] = G Q. In such a situation 
all stationary states satisfy well defined transformation 
laws, associated with an irreducible representation (ir- 
rep) r of dimension dp, which also specifies the respec- 
tive level degeneracy. For dp > 1, a degeneracy index 
j = l,...,(ir, the so-called "partner function", labels a 
choice of orthogonal states within the same eigenspace. 
Expressed for passive transformations, the laws read 



dr 



(33) 



i=l 



where D^{g) is a set of representation matrices that char- 
acterize the transformation laws of the partner function 
basis, and can be chosen in a suitable way. 

Since all stationary states can be associated with a well 
defined symmetry (r,j), the Hamiltonian can be pre- 
diagonalized by finding an orthogonal decomposition of 
the Hilbert space in subspaces gathering only states with 
symmetry (F, j). As for the one-phonon/one-electron QD 



model, this symmetry decomposition writes 



(34) 



where the orthogonal projectors Vv.j on the subspaces 
spanned by all the states that satisfy the transformation 
laws (|33]) for a given symmetry (F,j) can be written as 



(35) 



The problem of finding the QDPs reduces to solv- 
ing H^^^ inside each relevant subspace TYp ^ individ- 
ually. To provide these subspaces with suitable bases, we 
require a symmetrized eigenstate basis relative to ^ 
i.e. each basis state satisfies the transformation (|33]) for 
its particular symmetry (F,j). Such bases necessarily 
exist, since H'^ obeys the same symmetry diS H. To 
start with we symmetrize the electron subsystem and the 
phonon subsystem separately, i.e. 



electron : { l'?"^ ^)} 
phonon : { |0), |q)} 



{|Fe, je,Q^e>} 

{|0), l^phJph.C^ph)} 



(36) 



ae is usually a sequential index with energy, whereas aph 
represents a continuous degeneracy index because of the 
assumption of LO-phonon monochromaticity. The ex- 
plicit transformations ([36|) can be more subtle than an- 
ticipated. An example will be developed in detail for the 
Csy symmetry group in section llV Al These symmetrized 
bases allow the construction of a symmetrized basis of 
the tensorial products using generalized Clebsch-Gordan 
coefficients C^fjf^^^ ij^ the sense of point groups). 



(37) 



(a) |r,j;re,ae;0) = = T,j^ = j,ae) 

= X] ^'jfejlt I Te , Je , tte ) ® I Tp/i , jph , ttph ) 

Je J Jph 



Here F and j refer to the overall symmetry and Fg and 
Tph satisfy F C Fg Tp^. The phonon vacuum is always 
symmetrical, Tph = Ai, and hence the overall representa- 
tion of a state with zero phonons will always be identified 
with the electron representation, F = Fg (eq. [37|) . 

A symmetrized expression of the relevant subspace 7Y* 
immediately results from equation (|3Q|) by replacing the 
electronic energy index r with the pair index (Fe,ae). 
Expressing Hq^^ and Vr^,ae,i ph terms of the sym- 
metrized product basis (|37|) directly leads us to a set of 
non-orthogonal basis vectors, each of which transforms 
according to ([33|) for a particular symmetry (F, j). Those 
vectors can be regrouped in different "natural bases" B^j 
associated with the different subspaces defined in 

(|34|) . The expression of those vectors can be further sim- 
plified using the selection rule for the Frohlich matrix 
elements, which results directly from the transformation 
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laws and the invariance of the Hamiltonian, 

{r,j;re,ae;rph,aph\W"*\r',f;r'^ = r,ae;0) = 

unless (r, j) = (r, /) 

For the remaining non- vanishing matrix elements, we 
shah use the notation 



(38) 



(r, j; Te, ae, Tph, J?"* |r, j; T, a',; O) 
Finally, the natural bases B^ j write 



' |r,i;r,ae;0^, 



E: 



(39) 



The sum goes over ah indices, but one assumes that the 
Clebsch-Gordan coefficients C vanish, when Tph does not 
satisfy a selection rule F C Fg Vph- These bases are 
mutually orthogonal, but the vectors in each individual 
basis remain non-orthogonal. 

The overall symmetry index of an arbitrary natural ba- 
sis state is always equal to the symmetry index of the in- 
volved pure electron state (zero-phonon state), see equa- 
tions (|38|39p . Hence, if an existing representation F is 
absent in the considered set of bound electrons, there are 
no T-like polaron states^ even though we necessarily have 
F-like phonon states! 

Like in the previous two sections the natural bases ([39|) 
provide a prediction of the polaronic spectrum. Fig. [H 
shows the particular case of two bound electronic levels, 
where the second level is twice degenerate by virtue of the 
underlying dot symmetry. In particular, we emphasize 
the appearance of degenerate polaron levels, which can 
be associated with both, a degenerate or non-degenerate 
electron level. 

The dimensionality of the different subspaces 
be derived from the number of natural basis states for a 
given symmetry (F,j), 



dimiHrj) = card{B^j) = n^(l + n) 



(40) 



is the number of distinct electronic energies with a 
given symmetry F (o^g = 1^ • • • 7^^)7 and n denotes the 
total number of distinct electronic energies (n = n^). 
The dimension (|4Q|) is independent of the partner func- 
tion j in agreement with the feature that those functions 
can be defined arbitrarily inside a given representation 
F. Since X](r j) dim{l-t^j) equals dim{H*)^ given in (|32|) . 
expression (|^Q|) is a consistent refinement of the full di- 
mension. 



D. Computational Aspects 

Finding the polaron spectrum of the one-electron/one- 
phonon model reduces to diagonalizing the Hamiltonian 
jjQDP i]2side the low dimensional subspace spanned by 



QD potential 
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FIG. 4: (Color online) Qualitative structure of the polaron 
spectrum in the case of two electronic levels, one of which is 
twice degenerate due to the dot symmetry. The levels corre- 
spond to two different representations F marked in red and 
blue/green. The latter has two dimensions, characterized by 
the function j. Each gray bar indicates 3 more polarons that 
would result from interactions with two-phonon states. 



the natural basis S* given in ([28|) . ([3T]) or ([39|) . depend- 
ing on the physical situation. Although the number of 
Frohlich matrix elements for the interaction with LO- 
phonons has been minimized by the subspace reduction, 
their prerequisite computation can be numerically inten- 
sive for the arbitrary 3D wavefunctions that one should 
consider in a general case (see section HVl where a single 
wavefunction is typically sampled on 10^ points). To al- 
leviate this issue we have developed an original adaptive^ 
irregular discretization of the reciprocal space for lattice 
modes, and shown that it was an efficient method, also 
applicable when working directly with a non- orthogonal 
basis. 

The numerical benefit of an irregular reciprocal space 
discretization relies on the fast variation of the Frohlich 
matrix elements in function of the normal mode wave 
vector q in certain well localized domains. Increasing 
the local point density only in those domains remarkably 
improves the numerical precision with a minor increase of 
the required computational resources. To generate a well 
adapted irregular q-space discretization, we start with a 
regular coarse mesh covering the first Brillouin zone of 
the underlying lattice. Then the various Frohlich ma- 
trix elements are evaluated for all wave vectors q of the 
given mesh. This requires a preliminary computation 
of the volumes associated with each mesh node, taken 
as the volume of the respective Wigner-Seitz cells (Ap- 
pendix [VirC]) . The nearest neighbors with the highest 
difference between their Frohlich elements are added a 
new node in between, which provides the initial mesh for 
the next iteration. This algorithm is repeated until the 
maximal difference between neighboring Frohlich matrix 
elements falls below a preset threshold. Fig. [5] shows the 
Wigner-Seitz cells generated with this technique in the 
case of the first Brillouin zone of a body-centered cubic 
lattice. 

After the generation of the irregular q-space discretiza- 
tion and the computation of the respective Frohlich ma- 
trix elements, the set of natural basis vectors is evaluated, 
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(c) 































FIG. 5: (Color online) Polygonal Wigner-Seitz cells of the 
irregular q-space discretization, (a)-(c) zoom into first Bril- 
louin zone of a BCC lattice, (d) A few selected Wigner-Seitz 
cells from a region similar to (c). 



allowing to diagonalize the relevant Hamiltonian H^^^ . 

Let us finally evaluate the numerical value of working 
with a non-orthogonal basis. The sole consequence is that 
the standard eigenvalue problem becomes a generalized 
eigenvalue problem^ that is an equation of the type 



(41) 



QDP 



where \a) are the natural basis states, and Saa' is the 
so-called "mass matrix". This trade-off is advantageous, 
since optimized packages for the generalized eigenvalue 
problem are widely available, and one gets rid of an addi- 
tional basis change involving a Gram- Schmidt decompo- 
sition (often requiring enhanced precision for small scalar 
products). This is a positive numerical byproduct of the 
non-orthogonal theory. 

With the set of tools presented above, a spectral preci- 
sion down to O.OlmeV for typical QDPs can be reached in 
characteristic computation times of a few minutes using 
a present-day standard processor (3 GHz, 32 bit). The 
most computer intensive part is the preliminary evalua- 
tion of the Frohlich matrix elements. 



IV. APPLICATION TO PYRAMIDAL QDS 

A. Symmetrical Model and Non-Orthogonal Basis 

In this section, we apply the minimal model for the 
non-orthogonal theory (section to a realistic pyrami- 
dal GaAs/AlGaAs QD*2i This dot is part of a complex 



heterostructure represented by the geometrical model 
shown in Fig. [6^1.^^ We will take full advantage of the un- 
derlying Ca^ -symmetry group, which exhibits only three 
irreps Ai, A2 and E. The latter is two-dimensional 
and a possible basis results from symmetrizing £^-states 
with respect to the symmetry plane <Ji (spanned by the 
[111] and [112] crystalline directions in GaAs/AlGaAs). 
Thereby the partner function j = ± is identified with the 
parity index relative to ai. In graphical representations 
we shall consistently apply the color scheme: Ai (red), 
A2 (yellow), E-\- (blue), E— (green). QDPs will be com- 
puted using the symmetrized natural basis introduced in 
section IIII CI This basis will be derived analytically in 
three stages: (1) individual symmetrization of electronic 
and phononic eigenstates of (2) construction of a 
symmetrized product basis using Clebsch-Gordan coef- 
ficients, and (3) derivation of the symmetrized natural 
bases for the relevant subspaces j • 

First, we shall find symmetrized electron and phonon 
bases. As for the bound electron, all eigenstates of are 
automatically symmetrized and hence the task reduces to 
finding these eigenstates. This was recently achieved by 
Michelini et al.?^ using an effective mass model. For a 
dot height h = lOnm, there are two Ai-like levels (non- 
degenerate) and one E-like level (twice degenerate) as 
shown in Fig. [6)3. In the standard notation of section 
IIII CI i. e. iFg je-,OLe)^ those states write 



{|^i,l),|^±>,|^,2)} 



electron basis 



(42) 



where the index jg has been omitted in the case of the 
one-dimensional Ai -representations and the index ag has 
been omitted for the unique £^-level. We note that there 
are no A2-like electron states at low energy, which imme- 
diately predicts that there will be no 742-like QDPs (sec- 
tion [IIIC]). For the phonons (taken as bulk phonons) the 
symmetrization is inasmuch different as the eigenstates 
of H^^ such as plane waves ^q(x), are not automatically 
symmetrized. This feature relies on the monochromatic- 
ity assumption rendering all normal modes degenerate. 
A symmetrized eigenstate basis is properly derived in 
Appendix IVII Dl The resulting basis states superpose 
six (or four) plane waves, such that the directions of the 
different wave vectors are mutually related by symmetry 
operations (see Fig. [6t). We shall label such states with 
the respective vector q belonging to the subset which 
constitutes a sixth of the reciprocal space. In the case 
of £^-like superpositions there are two orthogonal states 
associated with the same vector q. They will be distin- 
guished through the additional index x = 1, 2 (discussion 
in Appendix IVII Dp . Finally the phonon basis writes 



0),|Ai,q),|A2,q), 1 



phonon basis 



(43) 



where |0) is the phonon vacuum state (OmeV), whilst all 
other states are one-phonon states (35.9meV). 

Second, we construct a symmetrized product basis 
from (US]) and (gS]) according to Eq.([37D. The explicit 
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FIG. 6: (Color) (a) Numerical modeP^ of realistic het- 
erostructure with pyramidal QD, vertical quantum wire 
(VQWR), vertical quantum wells (VQWs) and lateral quan- 
tum wells (LQWs). (b) Isosurfaces of envelope functions of 
the stationary single electron states, (c) Particular plane wave 
superpositions of symmetrized one-phonon basis. The length 
of the arrows represents the relative amplitude and dashed 
arrows have opposite phase. 



derivations given in Appendix IVII El yield states of the 
form 



ri;re 
ri;re 



=;0), 1 



product basis (44) 



where the first ket represents states with zero phonons 
and the latter states with one phonon. 

Third, we write the natural basis of the relevant 
Hilbert subspace of 7Y* according to the general the- 
ory (section fill Cp . This basis decomposes in symmetry- 
subspaces (8 dimensions), H%_^ (4 dimensions), 

^E- dimensions). Yet, in our particular case the most 
energetic natural basis states yield energies above the 
first two-phonon state. To remain consistent with the 
one-phonon assumption, we shall neglect those states. 
Thereby the dimensions reduce to 6 (H^J, 3 {Tl^_^) 
and 3 (7Y|;_). The respective natural bases result from 
the general expressions (|39|) and are given in Tab. [J and 
Tab.ini 



meV 


expressed as symmetrized product states 


subspace 


43.1 


|Ai;Ai,l;0) 




79.0 




Ai;Ai,l;Ai,q)(... 




Ai;Ai,l;0> 


108.1 








Ai;Ai,l;0) 


84.9 


|Ai;Ai,2;0) 




79.0 




Ai;Ai,l;Ai,q)(... 




Ai;^i,2;0) 


108.1 








Ai;Ai,2;0) 



TABLE I: Natural basis states of the subspace H\ . The bra 
(...| is the adjoint of the preceding ket. The two subspaces 
^Ai,i cind 'Hai,2 cire defined by the three basis vectors on 
their left. They constitute so-called "weakly coupled strong 
coupling regimes" , discussed in section IIV CI 



meV 


expressed as symmetrized product states 


72.2 


|^±;^;0) 


79.0 


Ea,A ^±;Ai,l;^,q,A)(... H^^' E±-E-0) 


108.1 


^r,„q,x ^; Tp,,, q, x) (..• |^±; ^; 0) 



TABLE II: Natural basis states of subspaces 7i^+ and TVe- 
The bra (...| is the adjoint of the preceding ket. 



B. Stationary States and Strong Coupling 

The problem of finding the stationary dot states, 
i.e. QDPs, consists in the eigenvalue problem 

i7^^^|rj,m) = er,.n|rj,m) (45) 

where m is a sequential energy index inside a particular 
symmetry (F, j). This eigenvalue equation was solved in- 
dividually inside each of the three decoupled subspaces 
^Ai' ^1;+ ^"e- using the enhanced matrix diago- 
nalization method outlined in section IIII Dl The three 
resulting spectra are given in Fig. [71 (red, green, blue). 
A geometrical representation of the corresponding po- 
laron states is shown in Fig. [HI where the closed surfaces 
are isosurfaces of the electronic and vibrational probabil- 
ity density functions. Those functions were obtained by 
computing the respective partial traces, 

Platticei'X.) = (x|Tre^ectron (|F j,m)(Fj,m|^ |x) (46) 
Pelectroni^) = (x| Tr/aUice ( | F j, m) (F j, m | ^ | x) (47) 

The two-dimensional representation E necessarily ex- 
hibits a spectrum consisting of twice degenerate levels, 
each of which is associated with one state in and 
one state in 



Each superposition c+|£^+,m) + 
C-\E—^m) is again a stationary state. 

Both the ground level and the first excited level yield 
negative energy shifts. This is consistent with the gen- 
eral feature that the ground level of each representa- 
tion is necessarily lowered with respect to correspond- 
ing free level. The numerical values of these shifts are 
Ae = —2.0 meV and = —2.5 meV. The same 
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FIG. 7: (Color) Spectrum of low energy quantum dot polarons 
(QDPs). Inside each symmetry subspace the states have been 
labeled with an energy index m, such that m = 1 corresponds 
to lowest energy state of a given symmetry. 



shifts computed with 2nd order perturbation theory are 
As = —7.0 meV and As = —17 A meV ^ respectively. 
This manifest large failure of a perturbative approach 
clearly confirms the existence of a strong coupling inside 
both irreps (Ai, E). 



C. Coupling Substructure 

For further characterization of the coupling regime it 
is interesting to consider the two state sets Si and 82^ 
defined as 



51 = {|Ai,l),|Ai,3),|^i,6)} 

52 = {|Ai,2),|Ai,4),|^i,5)} 



(48) 



The Ai -states of Fig. [8] have been ordered according to 
these sets. We demonstrated numerically that states in 
«Si are to a good approximation contained in the subspace 

1 ^ defined in Tab. [H Indeed, the norms of 

their projections on ^ exceed 95% of the full norms. 
With the same accuracy the states in ^2 are contained 
in ^Ai,2 ^ ^Ai' other words, the two subspaces 

I and 7Y^^ 2 appear reasonably decoupled^ although 
the whole subspace = 7i^^ ^ ^Ai 2 constitutes 

a strong coupling regime. Therefore the strong cou- 
pling regime must reside inside the two subspaces 7Y^^ ^ 
and 7Y^^ 2 individually, and they may be referred to as 
"weakly coupled strong coupling regimes" . The physical 
reason for this particular structure relies in the geome- 
try of the vibrational density function piatucei'^)- Fig. [8] 
shows that states in Si have a vibrational component, 
which is vertically centered in the dot, whereas the states 
in S2 have two centers of vibration splitting the isosur- 
face in two parts. Indeed the subspace i is spanned 
by two one-phonon states with vertically centered vibra- 
tional density and one zero-phonon state with centered 
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FIG. 8: (Color) Geometrical representation of QDPs. The 
two right columns show isosurfaces of the electronic and vi- 
brational probability density functions in direct space. 



electronic density. The resulting overlap leads to a strong 
interaction between electrons and phonons. The same 
conclusion applies to the subspace 2^ where the den- 
sity functions are vertically split in two parts. One the 
other hand, this picture reveals that the mutual overlap 
between 1~L\^ ^ and 2 considerably smaller. 

The concept of weakly coupled subspaces ^ and 
2 provides a direct tool for interpretation of the spec- 
trum in Fig. [3 In particular, the ground levels of each 
subspace, i.e. and |Ai,2), are necessarily lowered 

relative to the corresponding free levels. Analogically, 
the most excited levels of each subspaces, i.e. 5) and 
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1^1,6), are both raised. Their mutual sphtting remains 
very smah as they are to a good approximation uncou- 
pled. 

Finally, we emphasize that the novel concept of 
"weakly coupled strong coupling regimes", represented 
by the subspaces ^rj,i ^rj,2v • ^^^Y general and 
potentially applicable to all QDs. If the matrix element 
integral is close to zero due to the mutual orthogonal- 
ity of the electronic wave functions, these subspaces can 
be treated as decoupled in a good approximation. This 
idea is straightforward when working with the natural 
basis, and thus represents a further advantage of using 
non-orthogonal basis states. 

D. Entanglement measure and strong coupling, 
decoherence and relaxation 

An alternative characterization of the strong coupling 
is reflected in the entanglement of stationary dot states, 
i.e. QDPs. We have computed the entanglement between 
electronic and and phononic coordinates using the stan- 
dard measure introduced by Bennett et al..^^ For pure 
states, 

N 

Ent{\rj,m)) =-Y,hflogNhf (49) 

where the coefficients q are given by the diagonal 
Schmidt decomposition 

N 

|rj,m) = ^Ci\ei) (g) \phi) (50) 

i=l 

It is known that this form always exists for every par- 
ticular ket |r j, m), but the computation of the |rj, m)- 
dependent orthogonal vectors and in the 

electronic and phononic Hilbert spaces respectively, now 
requires a preliminary orthogonalization of the quantum 
dot phonon basis 5^^/|0) which can be elegantly per- 
formed via a Choleski decomposition of the A matrix. 
A subsequent singular value decomposition (SVD) of the 
coefficient matrix in the tensorial product basis will de- 
liver Ci . It is remarkable that the sum over i in ([5Q|) can be 
limited to the number of electron states because of the 
properties of the SVD. The entanglement measure (|49|) 
can vary between (non-entangled) and 1 (fully entan- 
gled), and is also often equivalently called the "entropy 
of mixing" of the two subsystems in state |rj, m). 

Tab. mil shows the entanglement of the QDPs pre- 
sented in Fig. [HI Weakly entangled states {Ent ^0.1) are 
nearly simple product states of electrons and phonons. 
In the present case, such a picture applies to the two 
QDPs 1^1, l) and |74i,3). Numerically, they consist 
to 99.5% of two natural basis states, which both in- 
volve the same electronic state |^i,l) (first two basis 
states in Tab. [Tj). All other QDPs are strongly entan- 
gled {Ent > 0.1) with no adequate perturbative pic- 
ture. Particularly strong entanglement is found in the 



states 1^1,5) and |yli,6), which involve peculiar Bell- 
state superpositions inside the ^^-representation of the 
form \E ^) ^\E ^) ^\E -) ^\E -). This of course 
suggests a general tendency to find particularly strongly 
entangled polarons in dots with symmetry related degen- 
eracies. 



State 


Entanglement 


Phonon Number 


Relax. Factor 




^1,1) 


0.023 


0.042 


0.001 




Ai,3} 


0.007 


0.963 


0.007 




Ai,6} 


0.523 


0.995 


0.520 




Ai,2) 


0.203 


0.920 


0.186 




Ai,4) 


0.222 


0.084 


0.019 




Ai,5) 


0.517 


0.996 


0.515 




E±,l) 


0.256 


0.136 


0.035 




E±,2) 


0.245 


0.902 


0.221 




E±,3) 


0.229 


0.962 


0.220 



TABLE III: Entanglement, 1-phonon probability and heuris- 
tic "relaxativity measure" of the Quantum Dot Polarons 
(QDPs) in the pyramidal dot with h = lOmn. 



Finally, we address a possible connection between en- 
tanglement and phonon-mediated decoherence and relax- 
ation. A simple model of such decoherence and relaxation 
would account for a weak bulk interaction between cou- 
pled LO-phonons and uncoupled LO-phonons or between 
coupled LO-phonons and LA-phonons. Such weak inter- 
actions maybe treated perturbatively and typically result 
in a finite lifetime for QDPs, which would otherwise be 
everlasting. One may expect that the lifetime depends 
on the entanglement between electronic and phononic 
coordinates, since entanglement indicates strong quan- 
tum correlations that could effectively translate phonon- 
phonon interactions to electron state hoppings. From 
this picture, we expect that the lifetime also scales with 
the weight of the phonon component. Hence, we heuris- 
tically propose a "relaxativity measure" for QDPs, de- 
fined as the product of the entanglement and the average 
phonon number 

Rel{\rj,m)) = {rj,m\Y,{^~^)x'xB{By\Tj,m) 

XX' 

xEnt{\Tj,m)) (51) 

In the present one-phonon model this measure varies be- 
tween (everlasting) and 1 (short coherence time, say 
~ Ips). At thermal equilibrium, the dot state is repre- 
sented by a density matrix exhibiting a high probability 
of states with a low relaxativity measure and vise versa. 
This measure does not really measure the relaxation since 
relaxation also implies other factors like resonances and 
population of final states, this is why we speak of "re- 
laxativity" . It has the status of a rough heuristic guess, 
since it is not the result of a proper relaxation model 
describing realistically phonon-phonon interactions and 
in particular neglects any dependence on the particular 
geometry. 
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E. Dot size variation V. INSIGHTS ON THE LOW ENERGY 

SCHEME IN QDS 



We shall now discuss the variation of the polaron spec- 
trum as a function of a varying dot height h. Fig- 
ure [9] shows the varying spectrum of (free levels) and 
the spectrum of H^^^ (quantum dot polarons) in the 
restricted energy band [40 meV, 120 meV] (quadrati- 
cally extrapolated from explicit computations of the dot 
heights 10 nm, 7.5 nm and 5 nm). To gain clarity and 
to remain consistent with the disappearance of certain 
free levels for smaller dots, we have restricted the graph 
to the two lowest levels of the two subspaces i and 
7Y|;^. The latter is of course degenerate with 7i^_ and 
the respective levels are twice degenerate. 




Dot height h [nm] 

FIG. 9: (Color online) Spectrum in function of the dot height 
h. (dashed lines) free electronic energies obtained by diago- 
nalizing , (solid lines) polaron energies obtained by diago- 
nalizing H^^^. 



There are three relevant free energies, the electronic 
ground level with symmetry Ai (red dashed line), the 
first excited level with symmetry E (blue dashed line), 
and the electronic ground level combined with one 
phonon (black dashed line). The latter two undergo a 
crossing in the vicinity of the dot height h = 2.5nm. 
By virtue of the resulting resonance, the two £^-like po- 
laron levels (blue solid lines) exhibit maximal energy 
shifts around h = 2.bnm (Ae ^ 5 meV) giving rise to 
a level anti-crossing. On the other hand, the decreasing 
resonance for increasing dot height, leads to a true cross- 
ing between the first excited E-\ike polaron level (upper 
blue solid line) and the first excited Ai-like polaron level 
(upper red solid line). Such a true crossing is consis- 
tent with the strict analytical decoupling stemming from 
group theoretical arguments {i.e. different irreps). 



We shall now expand the results to a very general 
class of QDs, including pyramidal, spherical, cubic or 
even cylindrical ones. For all these systems we uncover 
an analogous low energy spectrum, clear connections be- 
tween polarons and free levels, symmetry properties and 
qualitative dot size dependencies. 

Explicitly, we consider all dots with a non-degenerate 
electronic ground level and a twice degenerate first elec- 
tronic excitation. These dots include the special but pre- 
dominant class of dots with Cnv symmetry with n > 3. 
Qualitatively, they yield a low-energy polaron spectrum 
consisting of two shifted electron levels and a splitted 
electron+phonon level, see Fig. \W\ Group theory reveals 
three independent substructures (red, green, blue), where 
two are mutually degenerate (green/blue). This struc- 
ture can be derived from the natural basis ([39|) . or may 
be obtained from the spectrum studied above (Fig. [7j) by 
suppressing all QDPs with higher energies or associated 
with the third electron level. 

The fine dashed lines in Fig. [10] link each QDP level 
with the free level, from which it would arise, if one could 
continuously turn on the Frohlich interaction. These con- 
nections are important for understanding the QDP spec- 
trum, as levels within the same representation (here the 
levels with the same color) generally repel each other 
under the interaction. The relative position of the first 
excited polaron in the invariant representation (here 
|Ai,3)) depends on whether the free electronic energy 
spacing is larger or smaller than the constant phonon 
energy e^o, see Fig. [TOk and Fig. [TOb . In case (a), 
the state |Ai,3) can fall between the second electron 
level (green/blue) and the electron+phonon level (dashed 
level). In (b), the same state lies always above the free 
electron+phonon level. One can generally pass from sit- 
uation (a) to (b) by a dot size increase rendering the 
electronic energy spacing smaller than e^o- 

Further, the links between free levels and QDP levels 
allow to predict the variation of QDP levels with varying 
dot size. In general, the shifts become larger as one ap- 
proaches the resonance, which is the transition between 
case (a) and (b). The varying spacing between the up- 
per two free levels leads to a changing shift of the two 
degenerate QDP levels (green/blue). These changes are 
represented by the vertical arrows for increasing dot size. 
The dependance gets reversed when passing from case 
(a) to (b), due to the anti-crossing at the resonance. The 
shift of the two symmetrical polaron levels (red) is dot 
size independent, because of their symmetry-decoupling 
from the moving free level (green/blue) and the constant 
phonon energy clo- 

These findings are very generic. For example they 
agree with the results of Verzelen et al!^^ for the case 
of cylindrical dots (height/radius=12/18). Case (a) is 
obtained for radii < 13 nm, whilst case (b) correspond 
to radii > 13 nm. Following our discussion, it is straight- 
forward to understand that in their case 1 1 ± ) lies below 
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|5'1>, that |2 ± > lies above |P ± 0> and that \Sl) can 
only lie below |P ± O) for radius < 13 nm. We also see 
that the shifts of \S0) and \Sl) do not depend on the 
dot size due to symmetry decoupling and the constant 
spacing between |S'0) and \Sl). 



(a) A 

^ |E±,2> 

e-(2x) 

e- + ph |Ai.3> 


(b) 

— |E±.2> 

e- + ph |Ai.3> 




^^^^ — |E±,1> 

Slo (const) 




e- (2x) 

eLO (const) 1 |E±,1> 


free levels QDPs 


free levels QDPs 



FIG. 10: (Color online) Fundamental level structure for QDs 
with Cn7;-like symmetry with n > 3. (a) electronic en- 
ergy spacing exceeds phonon energy e^o, typically smaller 
dots, (b) electronic energy spacing smaller than e^o, typ- 
ically larger dots. The arrows show the level changes with 
increasing dot size. 



VI. SUMMARY 

In this work we uncovered the substantial advantage 
of the direct use of non-orthogonal creation and anni- 
hilation operators to treat polarons in quantum struc- 
tures. Starting from a general viewpoint, we fully re- 
formulated the polaron problem in terms of those opera- 
tors that naturally appear in the interaction Hamiltonian 
and generate the phonons relevant for individual transi- 
tions between electronic eigenstates. We also provided a 
complementary basis for all non-coupling phonons, which 
play a sensitive role in relaxation processes mediated by 
phonon-phonon interactions. Even though one might a 
priori be skeptic with the use of non-orthogonal objects, 
this approach proved mathematically elegant and fruit- 
ful for physical insights. In particular, we found a nested 
structure in the electron-phonon coupling, which allowed 
us to identify a non-trivial rule to truncate the Hilbert 
space in the case of a finite number of phonons. This fea- 
ture was consistently applied to a general QD structure 
in a one-electron/one-phonon model, and lead to a novel 
polaron basis, baptized the "natural basis". The latter 
constitutes an efficient tool for computation and detailed 
classification of quantum dot polarons (QDPs). Beyond 
the case of general quantum dots, we also investigated de- 
generate and symmetrical quantum dots using the appro- 
priate mathematical instruments, namely group theory. 
This revealed additional simplifications, degeneracies and 
subclasses of QDPs. 

As a realistic application we computed the low-energy 
QDPs of recently manufactured pyramidal QDs with 



-symmetry. To this end an adaptive irregular dis- 
cretization of the lattice mode space was developed, 
which we used to compute the Frohlich matrix elements. 
The generalized eigenvalue problem stemming from the 
direct use of non-orthogonal basis vectors was directly 
fed into efficient matrix diagonalization software. In this 
way, the requirement for computational resources was re- 
markably decreased. The numerical results explicitly re- 
vealed the spectral structure predicted from the natural 
basis. 3D-visualizations of the stationary polaronic dot 
states gave insight in the localization of both electronic 
and phononic components and showed the different sym- 
metry properties. Dot size dependent spectral investiga- 
tions uncovered level crossings and anti-crossings, which 
were consistent with the corresponding symmetry proper- 
ties. Further, we could prove the existence of strong cou- 
pling regimes for each symmetry representation through 
explicit comparison with second order perturbation the- 
ory. Yet, there was undoubtable numerical evidence for 
the presence of very weakly coupled subspaces within 
the strong coupling regimes. This led us to the con- 
cept of "weakly coupled strong coupling regimes" . Using 
the natural basis such subspaces could be understood in 
terms of specifically different overlaps between electronic 
wave functions and non-orthogonal vibrational modes. 
We used Bennett's entanglement measure to quantify the 
coupling between electronic and phononic coordinates — 
an idea that finally lead us to a heuristic "relaxativity 
measure" . In the end, we discussed the low-energy spec- 
trum of an important class of symmetric QDs (including 
spherical, cubic and cylindrical dots), and showed quali- 
tative predictions of the level structure and dot size de- 
pendance, valid as much for the general case as for our 
specific Csy pyramidal QD. 

We thank Pawel Machnikowski and our referees for 
their thorough and inspiring suggestions. Further, we 
acknowledge partial financial support from the Swiss NF 
project No. 200020-109523. 



VII. APPENDIX 
A. Derivation of the coefficients Cqq/ 
In section Hi B( the coefficients Cqq/ were defined as 

4 = E^qq'V^ (52) 

such that 5t|o) = (1 - V)bl^\0) with V being the or- 
thogonal projector onto vect{Bl^\o)}. By substitution 
we find 

Ecqq,6q,t|0) = (t_p)5^|0) (53) 

The projector V satisfies 

VBi\0) = Bi\0) (54) 
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and can be decomposed as 

V = Y,pxx'BI\0){0\B, 



(55) 



with complex coefficients pw. Substituting ([55]) in 
(f54j) and using the commutation relations (|6]) yields the 
unique solution p^A' = (A~^)aa'- Hence, 



r = Y,{A-')xx'Bi\0){0\B, 



XX' 



qq' AA' 



(56) 



Substituting ([56]) into ([53]) finally exhibits the unique 
solution 



Cqq' - (^qq' - ^)AA'^Aq'^^'q 



(57) 



AA' 



B. Demonstration of equation (|22p 

We want to find the p-phonon part H^^^ of a subspace 
<S, defined as 



vect< e 



(58) 



In the following, we implicitly assume that goes 
over all states of 7Y?ff (or, equivalently, over a basis 



of^p-T)- In the expansion of the exponential, the func- 
tions 1, t, t^, t^, . . . are linearly independent. Thus, 

S = vect\^{H^^^f\Lp) VA: = 0,1,...} (59) 

We then replace H^^^ by H^^W_^^^W_I^\ where H'^^ is 
the phonon creating term of and WJ^^ is the phonon 
annihilating term (^W^^ = H]^^ + i^!_^^), 

S = vect\^ {H^ + H!^' + HT')^ I (/:)) V A: = 0, . . . } (60) 

As we are interested in the p-phonon subspace of 5, 
the terms {H^ + i^^^ + j^int^k significantly sim- 

plified by retaining only the operator products increasing 
the phonon number by one unit. These are the products, 
which contain exactly one more H"^* than H^*. Fur- 
ther, we want to respect the assumed truncation of the 
phonon Fock space to at most p phonons, that is impos- 
ing = V G n^^^ (see section |irC|). Ex- 
plicitly, we need to remove all products involving inter- 
mediate {p-\- l)-phonon states (e.g. H^J^^W_^^ which 
involves the state H]^^'^\lp)). Applying these rules, the 
terms {H^ + + reduce to 



k = l 
k = 2 
k = 3 

etc. 



Since these vectors are used to span a collective subset, 
we can merely clean the list by creating new superposi- 
tions. Explicitly, we walk down the list from k = 1,2,... 
and subtract all the parts that are manifestly covered 
by smaller k already. For k = 2, for example, we can 
subtract W^^H^\ip) from {H^H_ 

H^\lp) = H'_^^\lp') with \lp') 



W^^H^)\ip), since 



e K 



QSP 



I is already 

spanned by the vectors associated with k = 1. Hence 
the additional vectors from k = 2 can be reduced to 
H^W^^\ip). {vect{W^'\ip)} and vect{H^W^^\ip)} are 
not necessarily linearly independent, but together they 
certainly span the same subspace as all the vectors in 
the list associated with k = 1 and k = 2.) We can then 
apply the same procedure to A: = 3 and find that all terms 
but H'^_^^\lp) are manifestly spanned by the vectors 
oi k = 1 and k = 2. One quickly realizes that proceed- 
ing in the same way, subsequently produces all the terms 
H^^W^*\ip), H^^W^*\ip), etc. Hence, 



n^^^ = vect\^{H'')^H!^'\^) V A: = 0, 1, ... } (61) 

Using again the property that 1, t, t^, t^, . . . are linearly 
independent functions of we finally find 



int I 



)Vi} 



(62) 



which concludes the demonstration. 



C. Expression of i/*^* for an irregular q-space 
discretization 

In the Frohlich matrix elements ([2]), the quantization 
volume V (direct space) dictates the underlying q-space 
discretization, such that each q occupies a volume of 
ft = {27r)^/V. This can be seen by taking V as a cubic 
volume with periodic boundary conditions, for which the 
Frohlich interactions was originally derived. If we use an 
irregular space discretization with varying cell sizes l^(q), 
the constant quantization volume V must consequently 
be replace by a function. 



V ^ V{q) 



(63) 



In the present case, ^^(q) was taken as the Wigner-Seitz 
volume around the point q in a given irregular reciprocal 
space discretization. 



C37;-Symmetrized Phonon Basis 

We consider the symmetry group Csv with its six group 
elements g = I (identity), g = C^, (positive and neg- 
ative 27r/3-rotation), g = (11,(72,(73 (plane symmetries). 
If |q) denotes the one-phonon state associated with the 
plane wave mode ^q(x), symmetrized one-phonon states 
|r,j, q) can be obtained by 

|r,i,q) = aPr,i|q)= ^ c^i(5)|7^(5)q) (64) 
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where 7Z{g) is the symmetry operation associated with 
the group element g. Vrj is the projector on the sub- 
space associated with the irrep F and the partner func- 
tion j. a is a normahzation factor defined up to a phase 
factor by the normahzation relations 



(r,i,q|r,j,q) = i vr,i,q 



(65) 



By projecting the rimodes basis states |q) on the four 
subspaces associated with Ai, A2, E+ and E—, one ob- 
tains an overcomplete set of 'irimodes states, which nec- 
essarily obeys irimodes relations of linear dependence. 
Those relations can be identified with the symmetry 
transformation relations, 



Ai:\Ai,n{g)q) = \Ai,q) \/g € Cs^ 



A2 ■■ \ A2,n{g)q} 



1^2, q) g = I,Ct,C3 
1^2, q) g = (Ji, (72,(73 



^. |E±,7^(<Jl)q) = ±|S±,q) 

■ \E±, q) + \E±, niC+)q) + \E±, niC^)q) 



(66) 



(67) 



= 
(68) 



q 

8 /A 


^1 







FIG. 11: Partition of the set of available wave vectors q. The 
whole set can be obtained by applying symmetry operations 
of the Csv group on the subset A. 



For the irrep Ai the five non-trivial relations of 
for a given q allow to restrict the plane wave set {q} to 
the sixth marked A in Fig. [Til The set {|Ai, q), q G ^} 
is orthonormal. An analog reasoning applies to the irrep 
A2 based on the five non-trivial relations of (|67|) . For the 
irrep E", (|68|) yields two non- trivial relations for a given 
vector q and a given partner function j = ±. Hence 
the set {q} may be restricted to a third of its elements, 
represented by ^US in Fig.[TTl Any two states j, q) 
and j, 7^(cr2)q) with q G w4 (and hence 7^(cr2)q G B) 
are non-orthogonal. In order to achieve orthogonality 
and to use one fixed vector set for all irreps, we introduce 
the states 

\E±, q, X = 1> = \E±, q) T ^2q) 
|E±,q,X = 2) = |E±,q)±|E±,a2q) 

where q G ^ and hence 7^(cr2)q G B. This definition 
completes the construction of the one phonon part of the 
phonon basis (|T3|) . The new index x = 1 , 2 permits the 



restriction of plane wave vectors to the sixth A and has 
the following physical interpretation: All ^-states with 
X = 1 involve one plane wave amplitude, whereas the 
£^-states with x = 2 mix two different amplitudes (see 
Fig.Efc). 

The canonical transformation relating the sym- 
metrized one-phonon basis to the plane wave basis {|q)} 
results from the definitions (|6T|) and (|69|) and the nor- 
malization relation (|65]) . We find. 



1^2, q) 

|^,+,q,i> 

|^,+,q,2) 

|^,-,q,i) 

\ |E,-,q,2)/ 



u 



^((T2)q) 
7^(<7l)q) / 



VqG^ (70) 



with the unitary transformation matrix 



u = 



1 

V6 

1 

V6 
1 
2 
1 

1 
2 

V 712 



1 

V6 

1_ 

V6 
_ 1 
2 

1 

^12 
1 
2 

1_ 

VT2 



1 

V6 
1 

V6 



1 



1 

Vs 



1 

Ve 
1 

Ve 


1 

Vs 



1 

Vs 



1 

Ve 
1 

Ve 
_ 1 

2 

1 

Vl2 
_ 1 
2 

1 

Vl2 



7-. ^ 

1_ 

1 
2 
1 

Vl2 
_ 1 
2 

"712 / 



(71) 



E. CsT^-Symmetrized Tensor Product Basis 

Based on the symmetrized electron basis (|^2]) and the 
symmetrized phonon basis (|43|) . we shall construct sym- 
metrized product states. The zero-phonon state |0) be- 
longing to the irrep Ai , the symmetrized product states 
involving zero phonons readily write. 



|Ai;Ai,ae;0) = \Ai,ae)^ |0> 
\E±;E;0) = \E±)^\0) 



(72) 



where n = 1 , 2 is the electronic energy index inside 
Ai. Semicolons separate intrinsic polaron, electron and 
phonon indices. As for the symmetrized product states 
involving one phonon, one uses Clebsch-Gordan coeffi- 
cients, 



\Ai;Ai,ae;Ai,q) 
\Ai;E;E,q,x) 

1^2; Ai,ae; A2,q) 
|^2;^;^,q,x) 

\E±;Ai,ae;E,q,x) 
\E±;E;Ai,ci) 
\E±;E;A2,q) 
\E±;E;E,ii,x} 



\Ai,ae) \Ai,q) 
■j^(\E+)®\E+,ci,x) 



+ \E-)<^\E-,q,x} 

\Ai,ae) ® 1^2, q) 

■j^(\E+)®\E-,ci,x) 

-\E-)®\E+,q,x)) 

\Ai,ae}®\E±,q,x) 
\E±)(g) \Ai,q} 

|^=f)® |^2,q) 

^{\E+)<E>\E±,ci,x) 

t\e-)®\Et,ci,x)) 
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where x = 1^2 is the additional phonon index used for 
^-hke one-phonon states. 
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